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' We present extensive tests of the Fast Action Method (FAM) for recovering the 

^ . past orbits of mass tracers in an expanding universe from their redshift-space coor- 

00 ' dinates at the present epoch. The tests focus on the reconstruction of present-day 

, pecuhar velocities using mock catalogs extracted from high resolution iV-body simu- 

' lations. The method allows for a self-consistent treatment of redshift-space distortions 

1*"^ . by direct minimization of a modified action for a cosmological gravitating system. 

' When applied to ideal, volume limited catalogs, FAM recovers unbiased peculiar ve- 

. locities with a 1-D, la error of ^ 220 Kms^^, if velocities are smoothed on a scale of 

' 5 h^^Mpc. Alternatively, when no smoothing is applied, FAM predicts nearly unbiased 

^ , velocities for objects residing outside the highest density regions. In this second case 

the la error decreases to a level of ^ 150 Kms~^. The correlation properties of the 
peculiar velocity fields are also correctly recovered on scales larger than 5 h~^Mpc. 
Similar results are obtained when FAM is applied to fiux limited catalogs mimicking 
53 i the IRAS PSCz survey. In this case FAM reconstructs peculiar velocities with similar 

^ ' intrinsic random errors, while velocity-velocity correlation properties are well repro- 

• i-H . duced beyond scales of ~ 8 h~^Mpc. We also show that FAM provides better velocity 

' predictions than other, competing methods based on linear theory or Zel'dovich ap- 

proximation. These results indicate that FAM can be successfully applied to presently 
available galaxy redshift surveys such as IRAS PSCz. 
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1 INTRODUCTION 

In the standard paradigm, the formation of cosmological structures is driven by gravitational amplification of tiny initial 
density fiuctuations (e.g. Peebles 1980). In addition to gravity, hydrodynamical processes can greatly influence the formation 
and evolution of galaxies, groups and clusters of galaxies. Hydrodynamical effects, however, play a minor role in shaping 
the observed distribution of galaxies on scales a few times larger than the size of galaxy clusters. Therefore, gravitational 
instability theory directly relates the present-day large scale structure to the initial density field and provides the frame- 
work within which the observations are analyzed and interpreted. Gravitational instability is a non-linear process. Analytic 
solutions exist only for configurations with special symmetry, and approximate tools are limited to moderate density contrasts. 
So, numerical methods are necessary for a full understanding of the observed large scale structure of the universe. There are 
two complementary numerical approaches. The first approach relies on TV-body techniques designed to solve an initial value 
problem in which the evolution of a self-gravitating system of massive particles is determined by numerical integration of the 
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Newtonian differential equations. Combined with semi-analytic models of galaxy formation, A^'-body simulations have become 
an essential tool for comparing the predictions of cosmological models with the observed properties of galaxies. Because the 
exact initial conditions are unknown, comparisons between simulations and observations are mainly concerned with general 
statistical properties. The second approach aims at finding the past orbits of mass tracers (galaxies) from their observed 
present-day distribution, independent as much as possible of the nature of the dark matter. The orbits must be such that the 
initial spatial distribution is homogeneous. This approach can be very useful for direct comparisons between different types 
of observations of the large scale structure. Most common are the velocity-velocity (hereafter v-v) comparisons between the 
observed peculiar velocities of galaxies are and the velocity field inferred from the galaxy distribution in redshift surveys (Davis, 
Nusser & Willick 1996, Willick et al. 1997, Willick & Strauss 1998, Nusser et al. 2000, Branchini et al. 2001) Ahernatively, 
one can also perform a density-density comparison between the mass density field inferred from the peculiar velocities of 
galaxies via POTENT-like methods (Bertschinger & Dekel 1989) or Wiener Filtering techniques (Zaroubi et al. 1995) and the 
observed galaxy distribution (see Sigad et al. 1998 and references therein) . Both types of analysis yield the cosmological mass 
density parameter Qm, for an assumed biasing relation between the distribution of galaxies and the mass density field. Any 
systematic mismatch between the fields serves as an indication to the nature of galaxy formation and/or the origin of galaxy 
intrinsic scaling relations used to measure the distances, provided that errors in the calibration have been properly corrected 
for. This second approach also allows to perform back-in-time reconstructions of the density field on scales ~ 5 h~^Mpc both 
in real (e.g, Nusser & Dekel 1992, Gramman 1993, Croft & Gaztaiiaga 1998, Frisch et al. 2001) and redshift space (Narayanan 
& Weinberg 1998, Monaco & Efstathiou 1999). Some of these methods has been applied to recover the past orbits from all-sky 
galaxy redshift surveys (Monaco et al. 2000, Narayanan et al. 2001). 

Finding the orbits that satisfy initial homogeneity and match the present-day distribution of mass tracers is a boundary 
value problem. This problem naturally lends itself to an application of Hamilton's variational principle where the orbits of 
the objects are found by searching for stationary variations of the action subject to the boundary conditions. The use of the 
Principle of Least Action in a cosmological frame- work has been pioneered by Peebles (1989) and has long been restricted 
to small systems such as the Local Group (Peebles 1990, Peebles 1994, Schmoldt & Saha 1998, Sharpe et al. 2001) and the 
Local Supercluster (Shaya, Peebles & TuUy 1995). Early applications to large galaxy redshift surveys have been hampered 
by the computational cost of handling the relatively large number of objects. Subsequent numerical applications speeded 
up the method and allowed the reconstruction of the orbits of ~ 10'^ particles (Shaya, Peebles & TuUy 1995). However, 
it was only recently that the improvement of the minimization techniques and the use of efficient gravity solvers made it 
possible to deal with more than 10* objects (e.g. the Fast Action Method [hereafter FAM] by Nusser & Branchini 2000 and 
the Perturbative Least Action Method by Goldberg & Spergel 2000), comparable to the number of objects contained in the 
largest all-sky galaxy redshift catalogs presently available such as the Optical Redshift Survey (Santiago et al. 1995, 1996) 
and IRAS PSCz (hereafter PSCz, Saunders 1996, Saunders et al. 2000). The situation has also been improved thanks to 
the invention of self-consistent direct schemes for accounting for redshift distortions, arising from the systematic differences 
between the distribution of galaxies in real and redshift space (e.g. Kaiser 1987). Previous applications (e.g., Shaya, Peebles 
& TuUy 1995) relied on iterations. However, self-consistent treatments of the problem have been proposed by Schmoldt & 
Saha (1998), Nusser & Branchini (2000), Phleps (2000), Susperregi (2001) and Goldberg (2001). 

In this work we will implement and test one of the Numerical Action Methods on realistic large mock redshift catalogs to 
reconstruct galaxy orbits over a large region of the universe. More precisely we will extend the Fast Action Method of Nusser 
& Branchini (2000, hereafter NB) to redshift space and focus on the problem of predicting galaxy peculiar velocities from a 
flux-limited, all-sky redshift catalog resembling the PSCz one. 

The outline of the paper is as follows. In Section ^ we review FAM in real and redshift space The performance of FAM 
is evaluated in the case of the ideal spherical infall in Section ^ Then a suite of more demanding tests is performed using 
the mock catalogs of galaxies described in Section ^. FAM is then applied to ideal, volume-limited mock catalogs in Section ^ 
and to fiux-limited mock catalogs in Section ^. Particular emphasis is given to the ability of FAM in returning an unbiased 
estimate of galaxy velocities and their correlation properties. The main conclusions are discussed and summarized in Section 
1^ We will use the terms FAMz and FAMr to refer to implementations of FAM in redshift and real distance space, respectively. 
Statements referring to FAM are meant to apply to both FAMz and FAMr. For brevity we will use x-space and s-space to 
refer to real-space and redshift-space, respectively. 



2 THE FAST ACTION METHOD 

In this Section we briefly summarize the Fast Action Method. For a detailed description of the method we refer the reader 
to Sections 2 and 4.2 of NB. We follow the standard notation in which a{t) is the scale factor, H{t) = a/a is the Hubble 
function, Q,m = p/pc is the ratio of the background matter density of the universe, p, to the critical density, pc = /SivG. 
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We denote the comoving coordinate of a patch of matter by x, and the corresponding comoving velocity by v = dx/dt. Also, 
let D(t) be the linear density growing mode normalized to unity at the present epoch, and f{Om) = dlnU/dlna ~ Q,^ (e.g., 
Peebles 1980). In the following we will use D{t) as the time variable in the equations of motion, with B — dx/dD as the 
corresponding velocity. Expressed in term of the time variable D, the equations of motion are almost independent on Sim and 
the cosmological constant (Gramann 1993, Mancinelli & Yahil 1995, Nusser & Colberg 1998). 

The evolution of a cosmological self-gravitating, isolated system of A'^ equal mass particles in a volume V is governed by 
the following equations of motion, 

de, 3 1 3 1 

dD^2W'-2D^jHSUn)^^^''' 

where the subscript i = 1 ■ ■ ■ A is the particle index, and g is the peculiar gravitational force field per unit mass. If the particles 
are unbiased tracers of the underlying density field then, 

I \ 1 \ ^ X — Xi 1 , , 

g(x) = --—y -. - + -X, (2) 

' 4™ ^ \x-Xi\^ i ^ ' 

i 

where n = N/V is the mean number density of particles inside V . In this expression for g we have assumed that NM — pV 
where M is the mass of a particle. These equations of motion can be derived from the action. 



+ — 

6 



(3) 



under stationary first variations of the orbits that leave x fixed at the present epoch and satisfy the constraint D^^^O 
as D (Peebles 1989, NB). The second condition on the velocities guarantees homogeneity as _D ^ 0. FAMr expands the 
orbits in a time dependent base functions g„ (D) in the form, 



x,{D) = a:,,o + V q„{D)C^.n, (4) 



where Xi^ is the position of the particle i at the present epoch, and the vectors Ci,„ are the expansion coefficients with 
respect to which the action is varied, i.e., they satisfy dS/dd^n = 0. The base functions and their derivatives p„{D) = 
dq^JdD are linear combinations of (1 — D), (1 — D)^ ■ ■ ■ (1 — /))""i£ix that the particle positions at _D = 1 is fixed, and 

ensures initial homogeneity. To simplify the expression of the gradient of the action, dS/dCi^n, 
NB imposed following the orthonormality condition on the functions p„{D), 



[ dDD'^^pJD)p,^{D) ^ 5^^„, 
Jo 



(5) 



where is the Kronecker delta function. The expansion of the orbits given in has to be modified when the constraints 
are galaxy redshifts, rather than the true distances. The redshift coordinate of an object is 

s„^Hox„ + {v„-s„)s„, (6) 

where the subscript refers to quantities at the present time, and So is a unit vector directed along the line of sight to the 
galaxy. When the constraints are the particle positions in s-space, S- „, the appropriate expansion of the orbits can be written 
as. 



n \ n / 



(7) 



where pn,o ~ Pn{D = 1), and the symbol (") indicates the component of a vector in the direction of line of sight at the 
present epoch, ^ . The trivial dependence on Ho can be completely eliminated by working with HqX instead of x. With this 
expression for the orbits, the gradient of the action subject to the new boundary conditions is 

as/aa,„ = i,,„ - /opn.o {e\^ + if , (8) 
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where /o = /(f^^.o), 



and 



Jo 



dBi 3 0i 3 1 r^Tn 



dfi[ 3e[_3_^an II 



(9) 



(10) 



The equations (^) differ from the equations of motion by the boundary term, /oP„ ^^1' ^ ^1 • This term can be 
eliminated by adding a kinetic energy term to the action, as follows 



5 = s+i^E(<o) 



(11) 



Given the redshifts of the particles, the minimization of the action S leads to the equations of motions (eq. 
transformations of this type have been proposed by Schmoldt & Saha (1998) and Withing (2001). 



Similar 



Our strategy is to find orbits that are as close as possible to the Hubble fiow. Therefore, we search for the minimum of the 
action and do not look for stationary points which might describe oscillatory behavior of the orbits (Peebles 1990, 1994). To 
find the coefficients Ci,n that minimize the action, FAM uses the Conjugate Gradient Method (CGM) which is fast and easy to 
implement (Press et al. 1992). The gravitational force g and its potential are computed using the TREECODE gravity solver 
(Bouchet & Hernquist 1988). The time integration in the expression for the action is done using the Gaussian quadrature 
method with 10 points at the time abscissa (Press et al. 1992). The CGM requires an initial guess for Ci,n. In the standard 
FAM application we compute the initial guess using the linear theory relation between the velocity and mass distribution. The 
minimum of the action proved to be rather insensitive to the choice of initial guess for Ci^n, as we have checked by running 
FAM experiments with initial Ci,n both set to zero and to random numbers with appropriate variance. Besides the initial 
set of Ci,„, the other free parameters are the softening used by the gravity solver and the tolerance parameter that sets the 
convergence of the CGM method. 



3 THE SPHERICAL INFALL MODEL TEST 

The collapse of a spherical over-dense perturbation can be followed analytically into the nonlinear regime before the occurrence 
of shell crossing (e.g., Peebles 1980). Testing FAM with the spherical collapse model will allow an assessment of the "shot- 
noise" errors in the recovered velocities. These errors arise from the discrete sampling of the density field and can be quantified 
by applying FAM in two different ways. In the first, we treat each particle as a uniform spherical shell. So the gravitational 
field on a particle is radial and only due to particles lying within its radius. In this way particles move along radial orbits and 
shot-noise errors are minimized. In the second, each particle is treated as point mass and the gravitational field is computed 
from the TREECODE as in any other application of FAM. The comparison between the peculiar velocities in these two 
different ways of applying FAM is an important test for assessing the performance of FAM under general conditions. 

Consider a spherical density fiuctuation described by the following radial density profile at the initial redshift Zin = 50: 



5in{r) = 5q 



(12) 



where ~ 40(1 + Zi^)^^ Mpc, Sq = (1 + 2i„)^^. The corresponding initial velocity is obtained using the velocity-density 
relation of the linear growing mode. We assume a fiat universe with = 1 and = 0. The analytic solution returns 
the density and peculiar velocity fields at a generic time. The evolved density field is then randomly sampled with 20,000 
particles within 80 Mpc. Each particle is then assigned a peculiar velocity according to the analytic solution and its redshift 
coordinate is computed accordingly. The distribution of these particles in x-space and s-space are then fed into FAMr and 
FAMz, respectively. We have used Umax = 6 base functions to expand the orbit of each particle, a softening of e = 0.5 Mpc, 
and a tolerance parameter tol = 10~* to determine the convergence of the CGM solution. 

In the first FAM experiments, particles move in radial orbits under the action of the radial gravity force. The final radial 
velocity profiles in the FAMr and FAMz experiments are shown, respectively, in the top-left and top-right panels of fig. ^ 
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Figure 1. Radial velocity profiles from FAMr (left) and FAMz (right) experiments. In all panels the continuous line shows the analytic 
solution for the spherical infall case. Filled dots represent the average FAM radial velocities measured in radial shells of 5 Mpc. The filled 
squares show the final overdensity profile, 5fi„(r) X 100. The numerical solution shown in the top panels is obtained by imposing radial 
FAM orbits. This constraint is removed in the solutions shown in the two bottom plots, where the scattered points represent FAM radial 
velocities measured at the position of a subset of particles. 



In both panels, the velocity profile in the analytic solution is represented by the continuous solid line. The filled dots show 
the average peculiar velocity computed within radial shells of 5 Mpc. The scatter around the mean is negligible. The filled 
squares represent the final overdensity profile, 5fi„(r) x 100, also in radial bins of 5 Mpc. Both FAMr and FAMz match the 
analytic solution very well at all radii. 

We have repeated the FAMz and FAMr experiments with the gravity force field computed by summing over the 3- 
dimensional discrete distribution of the point masses. This generates random shot-noise errors in the FAM recovered velocities. 
The radial velocity profiles are shown in the two bottom panels of fig. ^. The results of the FAMr and FAMz experiments are 
very similar. The average FAM solutions (filled dots) are still close to the true one. FAM appears to underestimate the true 
velocities in the high density regions. However, the effect is of little significance (~ 30 Kms~^) when compared to the 1-D 
scatter around the mean (~ 120 Kms~^) which, instead, does not seem to depend on the local density. 

The amplitude of the random errors is better appreciated in fig. |^ where we plot the Cartesian X-component of true 
vs. FAM velocities of ~ 1000 randomly selected particles. The parameters of the best linear fit and its 1-D, la scatter are 
shown in each panel. The diagonal line in each panel is plotted to guide the eye. The two upper panels show the results 
when the particles are forced to move along radial orbits. When this constraint is removed (lower panels), the shot noise 
error causes a 120 — 150 Kms^^ uncertainty in FAM reconstructions which is clearly visible both in x- (bottom-left) and 
s-space (bottom-right). The small systematic errors seen in fig. ^ are also visible in fig. |^ and cause the FAM velocities to be 
slightly underestimated, especially when the reconstruction is performed in x-space. Systematic errors can be quantified by 
the deviation of the slope of the best fitting line from unity. The main source for the systematic mismatch is the softening 
used in the TREECODE to compute the force field. As already pointed out by NB, the amplitude of FAM velocities decreases 
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Figure 2. Cartesian X-components of FAMr (left) and FAMz (right) velocities vs. true ones. Only velocities of 1000 randomly selected 
particles are shown. The upper panels show the case of purely radial orbits. The lower panels show the case of unconstrained orbits. The 
results of the least square fit and the 1-D, Icr scatter around the best fitting line are indicated in each panel. 

when increasing the softening parameter, along with the scatter around the best fitting Hne. This is true also in the spherical 
infaU experiment, as we have checked by running the same FAMr with a smaller softening parameter of e = 0.25 Mpc. As 
expected, using a smaller softening reduces systematic errors and increases the slope of the best fitting line from 0.92 to 0.95. 

We have also tested FAMz with spherical infall in a universe with Qm ~ 0.3 and f^A ~ 0.7. FAMz velocities are still free 
of systematic biases and the relative random errors are similar to those found in the previous tests in a flat, fim ~ 1 universe. 



4 THE MOCK iV-BODY CATALOGS 

In Sections ^ and ^ we will perform more realistic and demanding tests of FAM using a suite of mock catalogs extracted 
from two A''-body simulations performed by Cole et al. (1998). The initial conditions of the two simulations were generated, 
respectively, from two cosmological models of Cold Dark Matter: a flat ACDM model with flm = 0.3 and f^A = 0.7, and 
a flat rCDM universe with Qm = 1.0 and power spectrum shape parameter F = 0.25. Both simulations were ran with an 
AP'^M code loaded with 192^ particles in a box of side 345.6 h~^Mpc. In both simulations the initial amplitude of the density 
fluctuations was normalized to the observed abundance of galaxy clusters. This requires setting erg = 0.55^^"'* (Eke, Cole 
& Frenk 1996), where erg is the linear rms mass density fluctuation in top-hat spheres of radius 8 h~^Mpc (which hereafter 
we indicate with TH8). Before generating the mock catalogs the distribution of A''-body particles has been re-sampled by 
assigning a constant probability to each particle which is used to perform a Montecarlo rejection procedure from a Poisson 
distribution with mean equal to the required probability. The originaW-body particles particles are thus substituted by 0,1,2 
or on rare occasions an even higher number of particles that retain positions and velocities of the parent object. Particles at 
the same locations are usually collapsed into a single objects with higher mass. Occasionally, however, roundoff errors cause 
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Table 1. Main properties of the four sets of mock catalogs. Name: Label of the set. #Cat- Number of catalogs in each set (representative 
of ACDM and rCDM cosmological models), i^part- Number of particles in each catalog. Sampling: Monte Carlo sampling technique 
adopted. Kth- radius of the top-hat filter used to smooth the original simulation [l represents the average interparticle separation). 
Rmax- external radius of each sample. 



Name #cat #Part Sampling Rth ( h ^Mpc) Rmax { h ^Mpc) 



VL 5+5 20000 10000+10000 80 

VL5TH 5+5 20000 10000+10000 5 80 

FL 5+5 ~ 8000 PSCz Selection 80 

FL5TH 5+5 ~ 8000 PSCz Selection Max [5, 1] 80 



them to be retained as close pairs or triplets of particles. Furthermore, we have cooled down the particles' peculiar velocities 
since they are significantly larger than those of real galaxies. The new peculiar velocities v are obtained in two steps (Davis, 
Nusser & Willick 1996). First we perform a mass- weighted smoothing of the N-body velocities, VNbady, with a spherical top 
hat filter of radius 1.5 h~^Mpc and obtain a smoothed velocity field V1.5 Then the two vector fields are added linearly to 
give V = 0.7 X Vjvbody + 0.3 x vi.s. The new pairwise velocity dispersion measured for object with a relative separation of 
1.0 h~^Mpc is ~ 200 Kms^^, comparable to the value measured by Strauss, Ostriker & Cen (1998) for galaxies outside high 
density region in the Optical Redshift Survey. 

From these two simulations we have extracted 4 sets of mock catalogs. Each catalog lists particle positions, redshifts, and 
peculiar velocities in a spherical region of radius of 80 h~'^Mpc centered on a Local-Group look-alike particle (see Branchini 
et al. 1999 for details). Each set contains 5+5 independent mock catalogs, corresponding to the two simulations, so that in 
total we have 40 catalogs. The four sets are termed VL, VL5TH, FL, and FL5TH and their main properties are summarized 
in table 0. The sets VL and VL5TH are volume-limited catalogs, while FL and FL5TH are fiux-limited. 

A catalog in the VL set is generated by randomly extracting 20,000 particles in a sphere of radius 80 h~^Mpc, half of 
which sample the inner sphere of radius 40 h~^Mpc and the rest are contained in the shell between 40 and 80 h~^Mpc. 
Different masses have been assigned to internal and external particles to guarantee the same average mass density throughout 
the volume. A typical VL mock catalog is displayed in the top-left panel of fig. ^ which shows a slice of thickness 20 h~^Mpc 
cut through a VL mock catalog extracted from the rCDM simulation. 

To obtain the VL5TH catalogs we first smooth in x-space the particles' distribution in the A''-body simulation in a mass- 
weighted fashion with a spherical top-hat filter of radius 5 h~^Mpc and then we Poisson sample the smoothed density field 
with 20,000 particles, half of which inside 40 h~^Mpc. The bottom-left panel of fig. ^ shows the x-space particle distribution in 
one of the VL5TH catalogs. The particle distribution appears to be less clustered than in the unsmoothed case, as expected. 
The redshifts of mock galaxies are obtained by adding the radial component of the 5TH-smoothed A'^-body velocities to the 
particles' distance (expressed in Kms~^). 

The sets FL and FL5TH are designed to mimic the PSC« fiux-limited survey. A FL mock catalog typically contains 
~ 8000 particles within a sphere of 80 h~^Mpc, extracted from the A'^-body simulation according to the selection function of 
PSCz galaxies as determined by Branchini et al. (1999): 



<f>{r) 



1+^ 



-1.83 



if r < 6 h"^Mpc 
ff r > 6 h^^Mpc 



(13) 



where r is the distance from the center of the sphere in h~'^Mpc, ro — 6.05 h~^Mpc, and r* = 87 h~^Mpc. The mass 
assigned to each galaxy is equal to the inverse of the selection function at the galaxy position (t){x). An example of catalog in 
the FL set is shown in the top-right panel of fig. ^ The particle distribution is concentrated towards the center due to the 
cut in fiux. To obtain the FL5TH catalogs we follow a two step procedures. First we perform a mass-weighted smoothing of 
the particle distribution in the x-space within each of the FL mock catalogs using a top-hat filter with an adaptive radius 
of Rth = Max[5,Z] h~^Mpc, where I is the average inter-particle separation at the generic position. This smoothing radius 
compromises the need of avoiding nonlinear effects with that of minimizing shot noise errors. Then we sample the smoothed 
density with ~ 8000 particles with a radial distribution according to the selection function (eq. ^|). The particle distribution 
in one of the FL5TH mock catalogs is shown in the bottom-right panel of fig. ^ As for the VL5TH case, the redshifts of the 
mock galaxies in the FL5TH catalogs are obtained using the velocities smoothed with a top hat filter on the variable scale 
Rth- 
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Figure 3. Particle distribution in x-space in slices of thickness 20 h~^Mpc cut through mock catalogs in the tCDM A'^-body simulation. 
Each panel shows a VL (top- left), VL5TH (bottom-left), FL (top-right), and FL5TH (bottom-right) mock catalog selected from the 4 
sets. 



5 TESTS WITH VOLUME-LIMITED CATALOGS 

In this Section we apply FAM to the suite of VL and VL5TH catalogs. Unless otherwise stated, all FAM reconstructions 
were performed with a softening e = 0.5 h~^Mpc, a tolerance tol — 10"'', and Umax ~ 6 basis functions. Decreasing tol or 
increasing Umax makes little change to the final results. 



5.1 Unsmoothed Volume-Limited: VL 

In fig. ^ we show maps of the A'-body and FAMz peculiar velocities for one of the VL catalogs. The vectors represent the 
2D-projected peculiar velocities of particles in a slice of thickness 6 h~^Mpc cut through one of the VL-rCDM catalogs. 
The length of the vectors is drawn in units of 1 h~^Mpc = 50 Kms~^. The top-left panel shows FAMz velocities for all 
the particles in the slice. The A-body velocities of the same particles are shown in the middle-left panel, and the velocity 
residuals, Vjvtody — "Vfamz, in the the bottom-left. Points with the largest residuals are located in regions of high density. 
These regions are characterized by virial velocities which are not modeled correctly by FAM. In these regions, FAM typically 
predicts a large, coherent infall into the gravitational potential wells of density peaks. One example is the region centered 
around the point with (SGX, SGY) ~ (60, 0). The signature of virial motions is clearly visible in the A-body velocity map 
but is completely absent in FAM. Instead, FAM predicts coherent inward streaming velocities. 

FAM is meant to model gravitational dynamics in the mildly nonlinear regime. So when we move away from high density 
environments we expect an improvement in the agreement between the A-body and FAM velocities. To verify this hypothesis, 
we define Ss as the density contrast smoothed with a TH5 filter and we only compare velocities of particles in regions with 5s 
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Figure 4. Maps of 2D-projected peculiar velocities for points residing in a slice of thickness 6 h~^Mpc cut through a VL-tCDM catalog. 
The length of the vectors is drawn in units of 1 h~^Mpc = 50 Kms~^. The top row shows the FAMz-predicted velocities. Af-body 
velocities are shown in the middle row. The velocity residuals, Vjvbody — Vfaa/z, are displayed on the bottom. The maps shown in the 
panels to the left hand side refer to all the points in the slice while only the velocities of points with 85 < 3 and 85 < 1 are plotted in the 
central and right columns, respectively. 



less than a certain value. The central column of fig. |^ shows peculiar velocities of particles in regions with S5 < 3. Cutting at 
this density threshold excludes only '--^ 5 % of the particles, but considerably improves the agreement between the two fields 
(the "hot spot" located at (SGX, SGY) ~ (60,0) disappears from the map of the residuals). The agreement further improves 
when we exclude points in regions with (Js > 1 (~ 20% of the points) , as shown by the third row of panels in fig. ^ In this last 
case the velocity residuals (bottom-right) are further reduced and are not concentrated around regions of enhanced density. 
Very few large residuals still exist for points with large A'^-body velocities. 

Fig. |5| shows a more quantitative point by point v-v comparison. In each panel the Cartesian X-components of the A''-body 
velocities is plotted against the corresponding FAMz (top panels) and FAMr (bottom panels) peculiar velocities. Each panel 
shows the velocities of 500 particles randomly selected from the total of 50,000 lying within a distance of 40 h'^Mpc in all 
the 5 VL-rCDM catalogs. The left column shows the velocities of all particles. The panels in the middle and right columns 
show only the velocities of particles with (S5 < 3 and S5 < 1, respectively. 

To quantify the match between FAM and A''-body velocities we perform a linear regression of the X-Cartesian components 
of A'-body velocities on FAM assuming errors on the A'-body velocities only, using all the particles 50,000 in the 5 VL-rCDM 
catalogs. The average slope and 1-D scatter of the regression of the 5 catalogs in each set are shown in the corresponding 
panels of fig. ^ along with the rms values of the average slope and 1-D velocity scatter over the 5 catalogs. The reasoning 
behind performing the linear regression of A'-body velocities on FAM with errors on the A'-body velocities is twofold. First, 
as pointed out by Berlind, Narayanan and 'Weinberg (2000), A'-body velocities consist of a large-scale contribution predicted 
from the mass distribution plus an uncorrelated small-scale contribution that can be regarded as a thermal noise. This suggests 
putting all the errors on the TV-body velocities, while leaving the FAM-reconstructed velocities on the X-axis. Second, this 
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Figure 5. A'-body vs. FAM particle velocity scatter plots in x-space (bottom panels) and s-space (top panels). Only the X-Cartesian 
velocity compontent of 500 objects, randomly selected among the 50,000 contained within 40 h~^Mpc from the center of the 5 VL - 
tCDM catalogs are plotted. Particles in the left panels are selected from all density environments while only particles with < 3 and 
(55 < 1 are displayed in the central and right panels, respectively. In all the plots, we show the slope and the 1-D velocity scatter around 
the fit (which uses all the 50,000 objects available) while the straight line shows the expected correlation. 



choice mimics the real velocity-velocity comparisons in which the observed (A''-body) velocities are compared with theoretical 
(FAM) predictions and the regression is performed putting all the errors in the observed quantity. 

The top and bottom panels on the left hand side of fig. ^ show a clear correlation between A^-body and FAM velocities 
with small scatter around the best fitting line, both in x- and s-space, although in both cases FAM overestimates the amplitude 
of the velocities. A few out-liers are present in the plots. However, while FAMz (top-left) underestimates the amplitude of 
the velocities of the out-liers, FAMr (bottom-left) does the opposite. So the slopes of the linear fits in FAMr are smaller than 
the corresponding slopes in FAMz. When the comparison is restricted to particles with ^5 < 3 (central row) the number of 
out-liers decreases and the agreement with Ai'-body velocities improves significantly, indicating that discrepant velocities arise 
in high density regions as suggested by the velocity maps in fig. ^. This trend is confirmed by the scatter plots on the right 
column which only include points with ^5 < 1, well outside high density environments. In this last case the slopes of the FAMr 
and FAMz best fits are quite similar and even closer to unity. Furthermore, the out-liers have almost disappeared and the 
1-D velocity scatter around the fit decreased to a very small level of « 150 Kms~^, to be compared with a typical 1-D rms 
velocity dispersion of ~ 350 Kms~^ for the objects in the VL catalogs. We conclude that FAMz and FAMr perform similarly 
outside high density regions (i.e. for ~ 95% of the points). The results of the fits shown in the plots are summarized in the 
upper part of table |^. The column labeled by 4a, indicates the result of linear regressions in which out-liers in the v-v scatter 
plots are eliminated by iteratively discarding all points deviating by more than 4cr from the best fitting line. The 4(T-clipping 
procedure converges in a few iterations. The final outcome of this procedure is similar to that of imposing a density cut at 
5^ — 1 and returns a best fit closer to unity. 

So far only particles within 40 h~^Mpc were included in the v-v comparisons, although all particles out to 80 h~^Mpc 
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Table 2. FAMr and FAMz results for the VL (upper part) and FL (lower part) experiments for two different cosmological models (rCDM 
and ACDM). The average slope, b, and 1-D velocity scatter a (in Kms~^) are shown along with their rms scatter, measured over the 5 
independent experiments. The best fit values are listed for various overdensity thresholds <55 and for the case of 4cr-clipping. 









FAMr 








FAMz 






Model 




All 


55 < 3 


55 < 1 


4a- 


All 


<55 < 3 


^5 < 1 


4a 


VL-tCDM 


b 


0.58 ±0.04 


0.67 ±0.03 


0.78 ±0.03 


0.80 ± 0.03 


0.64 ± 0.03 


0.78 ±0.03 


0.82 ±0.03 


0.82 ±0.04 




a 


230 ± 24 


195 ± 11 


147 ±5 


131 ± 5 


253 ± 33 


214 ± 15 


161 ± 8 


141 ±8 


VL-ACDM 


b 


0.43 ±0.03 


0.71 ± 0.03 


0.84 ± 0.04 


0.88 ± 0.03 


0.68 ± 0.02 


0.82 ±0.03 


0.86 ±0.04 


0.91 ±0.06 




a 


287 ± 55 


161 ± 12 


124 ±7 


105 ± 10 


295 ± 59 


167 ± 12 


126 ± 6 


108 ±6 


FL-rCDM 


b 


0.67 ±0.03 


0.72 ±0.03 


0.78 ±0.02 


0.81 ±0.03 


0.74 ± 0.04 


0.80 ±0.04 


0.88 ±0.03 


0.92 ±0.05 




a 


227 ± 12 


201 ± 10 


171 ± 12 


152 ± 14 


242 ± 16 


212 ± 14 


170 ± 13 


145 ± 12 


FL-ACDM 


b 


0.51 ±0.02 


0.70 ± 0.03 


0.86 ± 0.03 


0.94 ± 0.03 


0.67 ±0.06 


0.86 ±0.04 


0.93 ±0.02 


0.96 ±0.03 




a 


223 ± 26 


168 ± 24 


135 ± 14 


108 ±9 


240 ± 60 


164 ± 24 


126 ± 15 


99± 11 



were used in FAM reconstruction. This is to guarantee dense sampling (the outer part of these catalogs is more sparsely 
sampled) and to minimize the eflFect of the external tidal field. Including all particles out to 80 h~^Mpc does not afTect the 
result appreciably, apart from a systematic ^ 5% increase in the slope of the best fitting line, which probably quantifies the 
bias introduced by having assumed a homogeneous mass distribution beyond 40 h~^Mpc in the former experiments. Similar 
considerations apply to the softening parameter. While, as shown in the NB analysis, increasing the softening parameter 
decreases the amplitude of FAM velocities, decreasing it from 0.5 to 0.25 h~^Mpc leaves the results of the v-v comparison 
unchanged. 

We have also checked the dependence of our results on the cosmological model by repeating the same analysis with 
the 5 VL catalogs extracted from the ACDM A'^-body experiment. The parameters of the linear fits to the v-v comparisons 
are also listed in table ^ and are similar to those found in the rCDM scenario. As for the rCDM case, FAMr velocities are 
systematically larger than the FAMz ones and the mismatch with A'^-body velocities decreases when moving away from high 
density regions. The effect is more dramatic in the ACDM universe, as it can be verified by comparing the results of the fits 
with no density cut with those of the 5$ = 3 cut. The reason for this behavior is that the ACDM model is characterized by 
a larger value of erg, i.e. by a larger number of regions with high overdensity where FAM fails to predict the correct peculiar 
velocities. 

The v-v scatter plots tell us little about the ability of FAM at recovering the correlation properties of the velocity 
field. These have been investigated using two different statistical tools. The first is the average relative pairwise velocity 
{V12) = {[v{x + ri2) — v(x)] ■ 612), where the averaging is over all pairs of objects at separation |t"i2|, and 612 is a unit vector 
along the the separation defined so that approaching pairs have (V12) > 0. The second statistics is the v-v correlation function 
projected along ri2, defined as {V1V2) = {v{x) ■ ei2v{x + ri2) ■ C12) (Gorski 1988). These statistics have been computed for 
each of the VL catalogs. The results are shown in fig. ^ for the rCDM case. Each "ribbon" represents the la uncertainty 
interval around the mean values of (V12) and (V1V2), averaged over the 5 reconstructions and plotted as a function of |i"i2|. 
The vertically-dashed and horizontally-dashed strips show the results of the FAMz and FAMr experiment. The dark strips 
show the "true" statistics computed from the A^'-body velocities. Similarly to fig. the two plots to the right have been 
obtained using all the particles within 75 h~'^Mpc while only particles with ^5 < 3 and ^5 < 1 have been used in the two other 
plots. 

The top-left panel shows that at separations smaller than ^ 1 h~^Mpc, both {Vi2)x and (Vi2)s are systematically larger 
than the A^-body result. This is the signature of what we have already spotted in the velocity maps, i.e. that in high density 
region (where most of the pairs with small relative separation reside) FAM predicts coherent infall ((V12) > 0) instead of 
disordered thermal motions ((V12) = 0). Also, the FAMr-predicted infall velocities are larger that the FAMz ones at separation 
smaller that few Megaparsecs. As we have already noticed in the linear fit to the v-v comparisons, this systematic discrepancy 
is caused by the smearing of density peaks in s-space which decrease the number of pairs having one member in over-dense 
regions and suppress the infall signal in the FAMz reconstructions. The same considerations apply to the plot on the bottom- 
left corner, which shows that the spurious infall makes {V1V2) of FAM velocities systematically smaller than the A'^-body one 
up to scales of ~ 3 h~^Mpc. These considerations are verified by the middle and right columns of fig. ^ Removing the high 
density regions does not change significantly the FAM average results but only the scatter. In contrast, the TV-body results 
change appreciably, especially at small separations. This indicates that most of the infall signal in the FAM velocities comes 
from close pairs (since the solution of first approach is the one preferentially found at the minimum of the action) which do 
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Figure 6. Relative pairwise velocity (V12) (upper panels) and velocity correlation function projected along the separation of the pair 
(yiV2) (lower panels), as a function of the separation r\2 averaged over the 5 VL-tCDM catalogs. The "ribbons" represent the la 
uncertainty intervals around the mean value. Vertically and horizontally-dashed ribbons show the results of the FAMz and FAMr, 
respectively. The dark strips refer to the A^-body experiment. The plots in the central row and those in the row on the right have been 
obtained by considering points with with is < 3 and is < 1 and within 75 h~^Mpc. All the particles available have been used in the 
plots on the left. 



not necessarily reside in high density environments. For the A*'-body experiment, however, the strong infall signature around 
5 h~^Mpc disappears when neglecting high density regions, meaning that most of the infall signal is contributed by pairs in 
which at least one member is located in a high density spot. 



5.2 Smoothed Volume-Limited: VL5TH 



The analysis in Section 5.1 revealed that when high density regions are excluded, FAM predicts an almost unbiased peculiar 
velocity field and reproduces its correlation properties on scales larger than a few Mpc. Instead of discarding regions where 
the velocity field is highly non-linear, we can apply FAM on a particle distribution obtained by sampling a smoothed version 
of the original particle distribution. Clearly, this smoothing operation is rather unrealistic, since real peculiar velocities are 
too sparse to be smoothed. Thus, the test performed in this Section is ideal. The aim is to minimize "thermal noise" and to 



5.1 



concentrate on the effects of sampling errors on FAM velocity prediction. Therefore, here we repeat the analysis of Section 
on the 5TH-smoothed, volume limited mock catalogs VL5TH. However, since the thermal noise has been eliminated by the 
smoothing procedure, in this section we assume that errors only affect FAM velocity prediction and perform linear regressions 
of FAM velocities on A'^-body, with errors on the Y-axis only. 

The plots in the central row of fig. |^ show the map of the 5TH-smoothed A''-body peculiar velocities interpolated at the 
particle positions in one of the mock catalogs, for points selected in slices similar to those shown in fig. ^ at the same three 
different density thresholds. The 5TH-smoothing operation has erased nonlinear signatures such as the virial motions in high 
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Figure 7. Peculiar velocities in a slice of thickness 6 h~^Mpc cut through one of the VL5TH-tCDM catalogs. Plots in the top row: 
FAMz velocities. Plots in the middle row: 5TH-smoothed A''-body velocities. Plots in the bottom row: velocity residuals Vjvbodj/ — 'Vfamz- 
Panels on the left: all the points in the slice are shown. Panels in the middle: points with ^5 < 3. Panels to the right: points with (S5 < 1. 



density regions. The A'^-body flow is dominated by large, coherent motions which are well reproduced by FAMz (upper row), 
as conflrmed by the small, incoherent residuals displayed in the plots on the bottom. As expected, none of the results in this 
Section change when applying different density-cuts. Therefore we will only show and discuss the results relative to the full 
sample. 

The scatter plots in fig. |^ show the v-v comparisons for the 5 VL5TH-rCDM experiments. The parameters of the linear 
fits are shown in each panel of fig. ^ and listed in table ^ The results of the FAMr experiments are very similar to the FAMz 
ones and do not depend on the cosmological scenario. In all experiments FAM velocities are slightly larger than the A'^-body 
ones. This discrepancy is due to the fact that the softening parameter in use (e = 0.5 h-^Mpc) is is much smaller than the 
smoothing radius. Increasing the softening to e = 3.0 h~^Mpc brings FAM and A'^-body velocities into agreement, as indicated 
in table ^. The 1-D velocity scatter around the fit is significantly smaller than the 1-D rms velocity dispersion in the VL5TH 
catalogs (Ri 240 Kms"'^). 

As shown in fig. ^, the analysis of the average relative pairwise velocity (upper panel) and the v-v correlation function 
(lower panel) demonstrates that the main differences between FAM and iV-body velocity fields are on scales smaller than 
~ 5 h~^Mpc. Indeed, at small relative separations (Vi2)s > (^12)2; and both are larger than the A'^-body ones, showing that 
FAM systematically over-predicts the relative infall velocity. This spurious infall pattern appears because of the small softening 
e = 0.5 h~^Mpc. As shown in fig. ^ varying the softening parameter has a large effect on the correlation properties of the 
FAM velocity field. With a smaller softening of e = 0.25 h~'^Mpc (horizontally dashed strips) the infall pattern becomes more 
prominent on small scales. Increasing the softening to e = 3.0 h~^Mpc (diagonally dashed strips) reduces spurious infall, and 
yields correlation properties which are very close to the A'-body ones over all scales. 
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VL5TH - tCDM-R < 40 Mpc/h 
All Particles 




Figure 8. FAM vs. A'^-body velocity scatter plots in x- (bottom panel) and s-space (top panel) for 500 objects randomly selected from 
the 5 VL5TH-tCDM catalogs. The straight line, showing the expected correlation, is plotted to guide the eye. The results of the linear 
fit are also displayed. 



Table 3. FAMr and FAMz results for the VL5TH experiments in the rCDM (upper part) and ACDM (lower part) cosmological models. 
The average slope, b, and 1-D velocity scatter, c, along with their rms scatter, are shown for FAM experiments performed with three 
different value of the softening parameter e = 0.25, 0.5, 3 h^^Mpc. 









FAMr 






FAMz 




Model 




e = 0.25 


€ = 0.5 


e = 3.0 


e = 0.25 


e = 0.5 


e = 3.0 


VL5TH-TCDM 


h 


1.07 ±0.08 


1.05 ±0.05 


1.08 ±0.10 


1.08 ±0.09 


1.05 ± 0.06 


1.02 ±0.09 






154 ± 22 


134 ±9 


112 ±8 


147 ± 13 


148 ± 10 


118 ±8 


VL5TH-ACDM 


b 


0.93 ±0.06 


1.03 ±0.05 


0.94 ±0.06 


0.93 ±0.10 


1.05 ± 0.06 


0.94 ±0.10 




a 


106 ± 21 


92 ±3 


97± 21 


115 ± 24 


98 ±6 


101 ± 24 



5.3 FAMz vs. Linear Theory and PIZA: Volume-Limited case 

In this Section we compare the FAMz velocities obtained in the VL experiments of Section to those predicted by linear 
theory and Zel'dovich approximation using the same VL-catalogs. Linear theory predictions have been obtained from the 
gravity force computed from the particle distribution in s-space at the final time, f Vi — -ffo/offo i- Note that no attempt 
is made here to account for s-space distortion effects which could be minimized using iterative procedures (e.g. Yahil et al. 
1991). Velocities in the Zel'dovich approximation have been obtained by running FAMz with Umax ~ 1 basis function which 
corresponds to straight line orbits Xi{D) — H^^S- „ + Dd.i — /oC"i,iS. These velocities should be similar to those that 
would be obtained by applying the PIZA method of Croft & Gaztanaga (1998) in s-space. 
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Figure 9. Relative pairwise velocity (V12) (upper panel) and velocity correlation function projected along the pair separation, (V1V2) 
(lower panel), computed from the 5 VL5TH-tCDM catalogs. The different strips represent la uncertainty intervals around the mean 
values. Vertically and horizontally-dashed strips refer to FAMz and FAMr experiments, respectively. The dark "ribbon" refers to the 
Af-body case. The correlations are computed using all the particles within 75 h~^Mpc. 



Fig. |ll| is similar to fig. g. The scatter plots in the upper panels show the comparisons between A''-body velocities and 
velocities predicted using linear theory (left), PIZA (middle) and FAMz (right) for 500 objects randomly selected from the 
50,000 particles in the 5 VL-rCDM catalogs. 

The scatter plot on the top left panel appears to be dominated by a large number of out-liers which cause linear theory 
predictions to overestimate peculiar velocities. This is not surprising and derives from sparse sampling, breakdown of linear 
theory in high density environments and from the fact that s-space distortions were not accounted for. To quantify the mismatch 
between A^-body and linear velocities we have performed a linear regression of A^-body velocities on model predictions using 
all particles in each of the 5 VL-rCDM catalogs and assuming errors on the Y-axis only. The average slope, 1-D scatter of 
the regression and their rms values are shown in the panel and listed in table ^ along with the results of the other velocity 
models. The slope of the best fitting line is very shallow and indicates that linear theory indeed severely overestimates peculiar 
velocities. The 1-D velocity scatter, however, is rather small compared to the 1-D A'-body velocity dispersion (~ 350 Kms^^). 
The situation improves considerably when using either the Zel'dovich approximation (top central) or FAMz (top right). In 
both cases the slope of the best fitting line increases dramatically, especially for FAMz. The scatter around the fit, instead, 
remains almost unchanged. 

Better results can be obtained by smoothing model velocities on small scales, where all the approximations break down. 
This is demonstrated by the scatter plots shown in the bottom panels of fig. ^ in which the unsmoothed A^-body velocities are 
compared with peculiar velocities predicted by linear theory (bottom left), PIZA (bottom central) and FAMz (bottom right) 
all of them smoothed with a Top Hat filter of radius 5 h~^Mpc. In all cases the match between true and predicted velocities 
improves considerably, especially for the case of linear theory. However, only FAMz returns unbiased velocity predictions. 
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Figure 10. The same functions (V12) (upper panel) and (Vi V2) (lower panel) displayed in fig. ^ are plotted here for the FAMz case and for 
three different values of the softening parameter. Horizontally, verticallv and diagonally dashed strips refer to e = 0.25, 0.5, 3.0 h~^Mpc, 
respectively. Note that the scales of the Y-axes are different from fig. bl 



Two sources of errors mainly affect velocity reconstruction methods: sparse sampling and the breakdown of model predic- 
tions in the highly nonlinear regime. Both uncertainties affect all velocity- velocity comparisons of fig. ^ (and those in fig. 3 of 
Nusser & Branchini (2000)), and therefore do not affect the relative comparison of the three reconstruction methods. However, 
these errors need to be accounted for to estimate absolute errors or when comparing our results with those of other tests that 
adopt similar velocity models but use different mock catalogs, as in Croft & Gaztafiaga (1998) and Berlind, Narayanan & 
Weinberg (2000). For example, the average inter-particle separation of 4.75 h~^Mpc in our mock VL catalogs is similar to the 
galaxy-galaxy separation in the inner part of redshift catalogs like PSCz, 2dF (CoUess et al. 2001) and SDSS (Friemann & 
Szalay 2000) but is larger than that of particles in mock catalogs used to test most velocity models. Furthermore, as explained 
in Section ^ our mock catalogs contain spurious pairs and triplets of nearby objects. Their presence artificially increases the 
amplitude of model velocities in high density environments and amplify the magnitude of the reconstruction errors. Both 
effects conspire in making our error estimates larger than those obtained from more ideal mock catalogs, like those used in 
the analyses of of Croft & Gaztafiaga (1998) and Berlind, Narayanan & Weinberg (2000). 

To check the dependence of our results on the cosmological model we have repeated the same analysis using the 5 VL- 
ACDM catalogs. The results of the fits to the v-v comparisons are listed in table ^. The performances of PIZA and FAMz in 
a ACDM universe are similar to those obtained in a rCDM one, indicating that the goodness of the velocity reconstruction 
depends little on the underlying cosmological model. The failure of linear theory, however, is more severe in the ACDM 
universe. This is a consequence of the fact that this cosmological model is characterized by a larger number of high density 
peaks in which linear theory breaks down. 

Smoothing FAMz velocities with a 5TH filter returns unbiased velocity predictions in both rCDM and ACDM scenarios. 
Increasing the smoothing length would have the effect of bring PIZA velocities (and, eventually, the linear ones) into agreement 
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Figure 11. A^-body velocities vs. peculiar velocities predicted by linear theory (left), PIZA (middle) and FAMz (right) for ~ 500 objects 
randomly selected from the 50,000 particles in the VL-rCDM catalogs. In the upper plots the comparison is performed using unsmoothed 
peculiar velocities. The predicted velocities in the lower panels have been smoothed with a 5TH filter and are compared with unsmoothed 
A'^-body velocities. The results of the fits are shown in all the plots. The straight lines at 45° are drawn to guide the eye. 



with the A'^-body ones. However, as pointed out by, Berhnd, Narayanan & Weinberg (2000), increasing the smoothing scale 
would also cause the FAMz velocity to be underestimated, as a consequence of the fact that errors in predicted velocities 
correlates with the predicted velocities themselves. We conclude that 5 h~^Mpc should be regarded as the optimal smoothing 
scale for FAMz velocities. 



6 TESTS WITH FLUX-LIMITED CATALOGS 

Real datasets typically consists of flux (or magnitude) limited catalogs. Here we investigate how well FAM performs when 
applied to the more realistic mock flux limited catalogs described in Section ^. We will be concerned with flux limited samples 
and will neglect other possible sources of systematic errors like unobserved regions of the sky (e.g. Zone of Avoidance), the 
occurrence of morphology and density biases and so on. 

In all experiments we have chosen to use the the selection function of PSCz galaxies (eq. since that catalog still 
represents the deepest all-sky redshift survey presently available. The mass of each mock galaxy is equal to the inverse of 
the PSCz selection function measured at the galaxy's x-space location. This may generate systematic errors when the FAM 
reconstruction is performed s-space (i.e. with FAMz) which are included in the error budget in the following two Sections. 
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Table 4. Af-body velocities vs. linear theory, PIZA and FAMz velocity predictions for the VL (upper part) and FL (lower part) 
experiments and for two different cosmological models (rCDM and ACDM). The average slope, b, and 1-D velocity scatter a (in 
Kms"'^) are shown along with their rms scatter, measured over all particles in the 5 independent experiments. The best fit values are 
listed for the case of unsmoothed predicted velocities (left hand side) and 5TH-smoothed predicted velocities (right hand side) . 









Unsmoothed 






5TH-smoothcd 




Model 




Linear 


PIZA 


FAMz 


Linear 


PIZA 


FAMz 


VL-tCDM 


b 


0.14 ±0.05 


0.58 ± 0.03 


0.64 ±0.03 


0.56 ±0.04 


0.84 ±0.02 


0.96 ± 0.02 




(7 


291 ± 31 


240 ± 30 


253 ± 33 


238 ± 30 


222 ± 14 


225 ± 27 


VL-ACDM 


b 


0.03 ±0.02 


0.51 ±0.07 


0.58 ±0.04 0.26 ±0.04 


0.81 ±0.05 


0.98 ± 0.06 






320 ± 50 


285 ± 40 


295 ± 60 


280 ± 57 


254 ± 35 


271 ± 56 


FL-rCDM 


b 


0.13 ±0.02 


0.52 ±0.07 


0.74 ±0.04 0.65 ±0.05 


0.93 ±0.04 


1.03 ± 0.04 




(J 


291 ± 23 


249 ± 25 


242 ± 16 


232 ± 18 


223 ± 11 


213 ±9 


FL-ACDM 


b 


0.04 ±0.01 


0.47 ±0.06 


0.57 ±0.06 


0.28 ±0.06 


0.88 ±0.04 


1.01 ± 0.04 




(J 


270 ± 60 


210 ±26 


232 ± 31 


223 ± 52 


192 ± 28 


209 ± 23 



6.1 Unsmoothed Flux-Limited: FL 

We have applied FAMz to the unsmoothed, flux limited catalogs labeled FL in table ^ and displayed in the upper-right panel 
of fig. ^. In this Section we present the results of v-v analyses analogous to those in Section 5.1. 

The v-v scatter plots for FL-rCDM catalogs are shown in fig. |l^, both for FAMz (upper panels) and FAMr (lower panels), 
for three different overdensity cuts. The scatter plots look similar to those of the volume limited tests (fig. The parameters 
of the linear regression of A'^-body velocities on FAM which are listed in the lower part of table ^ confirm the visual impression 
that FAM velocities are systematically larger than the TV-body ones. Moreover, FAMr produces systematically larger velocities 
than FAMz. The mismatch and between FAM and A'^-body velocities decreases outside high density environments. For < 1, 
FAMr and FAMz velocities are only 10—20% larger than the A-body ones and the mismatch further decreases when performing 
the 4(T clipping procedure described in Section 5.1, The analogy with the results found in is Section 5.1 shows that selection 



effects do not seem to affect FAM performances too much. The most striking feature is perhaps that the scatter around the fit 



is very similar to that found in the volume limited experiment of Section 5.1. FAM performances are similar in the rCDM and 
ACDM cosmological models despite of their very different values of ag, and greatly improve when high density environments 
are excluded from the analysis. 

Also the results of the velocity correlation analysis are quite similar to those of the unsmoothed- volume limited tests, as 
it is evident from the similarities between figs. ^ and However, selection effects are responsible for three main differences. 
First, the amplitude of the difference between FAM and A-body velocity correlations below ~ 1 h~^Mpc is twice larger than 
in the volume-limited case. Second, the agreement with the A-body results shifts to a somewhat larger scale. Third, the 
scatter around the mean correlation value increases. 

The results of this Section show that, when applied to mock PSCz fiux-limited catalogs, FAMz produces almost unbiased 
estimates of peculiar velocities outside high density regions. Moreover, the correlation properties of the FAM and A-body 
velocity fields agree on scales larger than 5 — 8 h~^Mpc. These results are remarkably close to those of the volume- limited 
tests. 



6.2 Smoothed Flux-Limited: FL5TH 

In the previous sub-section we have applied FAM to the particle distribution in the FL catalogs. Here we perform an ideal 



experiment analogous to that of Section 5.S and show the results of FAM applied to the FL5TH catalogs in which the objects 
sample a smoothed version of the particle distribution in the FL catalogs. A typical FL5TH mock catalog is displayed in the 
bottom-right panel of fig. |^. As for the case of the VL5TH tests discussed in Section. the results do not depend on the 
local density and therefore we only show the results relative to the full sample. 

Fig. shows the v-v comparison for the 5 FL5TH-rCDM experiments. As in Section. 5^ we have performed linear 



regressions of FAM velocities on A-body ones, where errors are associated to FAM velocities only. The results are displayed 
in fig. ^ and in table ^. They show that FAM peculiar velocities are systematically smaller than the A-body ones. Such 
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Figure 12. A'^-body vs. FAM velocity scatter plots in x- (bottom panel) and s-space (top panel) for 500 randomly selected objects from 
the 5 FL-tCDM catalogs. Panels on the left: all particles. Panels in the middle: particles with <55 < 3. Panels on the right: particles with 
5^ < 1. The straight line shows the expected correlation. The results of the fit are shown in all the panels. 



Table 5. FAMr and FAMz results for FL5TH experiments in the rCDM (upper part) and ACDM (lower part) cosmological models. 
The average slope, b, and 1-D velocity scatter, cr (in Kms~^), measured over the 5 independent experiments are shown along with their 
rms scatter. 



FAMr FAMz 


Model 




All All 


FL5TH-TCDM 


b 


0.93 ±0.10 0.91 ±0.11 




cr 


174 ± 17 288 ± 75 


FL5TH-ACDM 


b 


0.93 ±0.07 0.99 ±0.20 




cr 


142 ± 27 256 ± 134 



discrepancy, however, is below the la significance level. The random errors are larger than in all previous tests, especially for 
the FAMz method, probably indicating that the smoothing-resampling procedure enhances the amplitude of the shot noise 
errors. It is important to stress that, like in the VL5TH tests, the results do not change when applying different density cuts, 
which indicates that by applying the smoothing-resampling strategy there is no need of identifying high density environments 
and exclude them from the v-v analyses. FAMr results seem to depend little on the underlying cosmological model, as indicated 
by the results listed in the lower part of table ^ This is not true for FAMz, indicating that in a universe with a large value of 
as, characterized by the presence of many high density peaks, some of them enhanced by effect of the shot noise and smeared 
by s-space distortions, the smoothing-resampling procedure might not be sufficient in providing FAMz-velocity predictions 
which are uniform throughout the volume of the sample. 
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FL - tCDM - R<75 Mpc/h 




Figure 13. (V12) (upper panels), and (V1V2) (lower panels), as a function of the pair separation ri2 from 5 FL-rCDM experiments. 
Vertically-dashed strip: FAMz. Horizontally-dashed strip: FAMr. Dark strip: TV-body. Left panels: all the particles are considered. Central 
panels: particles with 5^ < 3. Right panels: particles with 5^ < 1. 

The effect of shot-noise is particularly dramatic for FAM velocity correlations. As seen in fig. the random uncertainties 
in the predicted velocity correlations become extremely large below the scale of ~ 2 h~^Mpc. On scales larger than 2 h~^Mpc 
the results look similar to those obtained from the flux limited tests, displayed in the plots on the left hand side of fig. |l2| 
In particular, the la uncertainty strip of the FAM-predicted (V12) statistics overlaps with the A'^-body result on scales larger 
than ~ 8 h~^Mpc. A similar behavior also characterizes the ly\V2) statistics. In this case, however, the agreement between 
FAM and A''-body is already found down to scales as small as ~ 3 h~^Mpc. 



6.3 FAMz vs. Linear Theory and PIZA: Flux-Limited case 



In this section we repeat the analysis performed in Section and compare the FAMz velocities discussed in Section 3.1 with 
those predicted by linear theory and PIZA in s-space, using the same FL catalogs. 

The results for the 5 FL-rCDM catalogs are displayed in fig. |l^ which is analogous to fig. |l^ and show the v-v comparison 
between unsmoothed A''-body velocities and velocities predicted by linear theory (top left), PIZA (top central) and FAMz 
(top right). Predicted velocities in the lower panels were smoothed using a Top Hat filter with an adaptive radius Rth = 
Max[5, l\ h~^Mpc, where I is the average particle-particle separation at the given redshift. The results of the linear regressions 
of TV-body velocities on the predicted ones are displayed in fig. |l6| and in the lower part of table ^. 

When no smoothing is applied to predicted velocities (upper panels), all methods overestimate the peculiar velocities. 



The situation is thus similar to that of the volume limited experiment performed in Section 5.3. This bias is particularly 
severe in the case of linear theory and decreases significantly by applying a top hat smoothing of radius Rth — 5 h~^Mpc to 
predicted velocities, as shown by the scatter plots in the bottom panels. In particular, FAMz smoothed velocity predictions 
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Figure 14. Same as fig. g but referring to the 5 FL5TH-rCDM mock catalogs. 

are bias-free and with a random error of of 220 Kms~^, nearly independent on the underlying cosmological model. These 
results demonstrate that, unlike the two other method considered here, applying FAMz to a PSCz-hke redshift catalog allow 
to predict a model for the cosmic velocity field which is unbiased, provided that FAMz velocities are smoothed on a scale of 
~ 5 h~^Mpc. 



7 SUMMARY AND CONCLUSIONS 

We presented extensive tests of the FAM method for reconstructing peculiar velocities from particle distribution in x- and 
s-space. Using the analytic spherical collapse model and a suite of realistic volume-limited and fiux-limited mock catalogs, we 
showed that FAM can successfully recover peculiar velocities from the particle distribution in x- and s-space. Errors in the 
predicted velocities are mainly due to shot noise and highly nonlinear motions that FAM fails to model correctly. Shot noise 
errors depend on the sparseness of galaxy sampling and thus on the particular dataset considered. To obtain a realistic error 
estimate we set the number density of galaxies in our mock catalogs similar to that in the inner part of the real PSCz, 2dF and 
SDSS catalogs. Errors due to nonlinear motions can be reduced by avoiding high density regions. Indeed, discarding objects 
in the highest density peaks brings systematic errors down to a level smaller than 20% and reduces the amplitude of random 
uncertainties from ~ 250 Kms~^ to ~ 150 Kms~^ in the volume limited experiments. When applied to flux-limited mock 
catalogs mimicking the distribution of galaxies in the PSCz redshift survey, FAM systematic errors outside dense regions are 
still below 20 % and the 1-D random velocity errors amount to 120 — 170 Kms^^. Moreover, FAM recovers the correlation 
properties of the underlying velocity field down to scales of ~ 5 h~'^Mpc in both flux and volume limited cases. In practice, it 
may be difficult to measure the local density from sparse datasets, in order to discard high density peaks. However, particles in 
these regions have very large velocity residuals, and appear as out-liers in the v-v scatter plots. Eliminating all points deviating 
by more than 4ct from the best fitting line improves the results of the fit and returns velocity estimates which are almost 
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F5TH - tCDM - R<75 Mpc/h 
All Particles 




(Mpc/h) 

Figure 15. Same as fig. ^ for the case of 5 FL5TH-rCDM experiments. 

unbiased. Another strategy to reduce the impact of nonlinear motions is that of smoothing the model velocity field. We have 
found that FAMz velocities smoothed with a Top Hat filter of radius 5 h~^Mpc are free of systematic biases. Clearly, there 
is nothing fundamental in this smoothing scale. To the contrary, the amount of smoothing required to minimize systematic 
errors depends on the strength of nonlinear motions, sampling rate and on the particular method used to compare model 
and observations. Other reconstruction methods like those based on linear theory and Zel'dovich approximation could also 
predict unbiased velocities, provided that an appropriate smoothing filter is applied. However, the smoothing scale required 
by FAMz to eliminate systematic errors is smaller than that required by linear theory or PIZA and thus allow to retain more 
information which lead to a more precise estimates of cosmological parameters such as f) — /b. 

Since FAM errors depend on several external factors, like the particular dataset considered, it is difficult to compare 
our results with those of other analyses. For example, it is somewhat surprising that errors in the Cartesian components of 
FAMz velocities predicted from our mock PSCz catalogs are a factor of two larger than those of the velocities obtained by 
applying the the iterative, quasi-linear model of Sigad et al. (1998) to the IRAS 1.2 Jy. catalog (see Willick et al. 1997). 
Indeed, the better treatment of nonlinear dynamics by FAMz and the denser sampling of the PSCz catalog should reduce 
these uncertainties rather than increase them. As discussed by Branchini et al. (2001), this apparent inconsistency derives 
from the differences in the mock catalogs used to calibrate the two reconstruction methods. Peculiar velocities in our mock 
catalogs, extracted from the AP^M simulations of Cole et al. (1998), are more nonlinear than those of the mock catalogs 
used Willick et al. (1997) that were extracted from PM A''-body simulations. The higher "temperature" of the velocity field 
in our mocks increases the uncertainties in the velocity reconstruction procedure. Errors in velocity predictions are further 
amplified by the spurious pairs and triplets of nearby particles contained in our mocks. The net result is that our estimated 
FAMz velocity errors are larger that those of Willick et al. (1997). Another example is given by the large systematic errors 
affecting our linear velocity predictions compared to FAM velocities, as shown by figures ^ and ^| This result suggests that 
estimates of /3 from v ~ v comparisons that use a linear model velocity field based on the PSCz catalog (e.g. Nusser et al. 
(2000) and Branchini et al. (2001)) are biased low and that this bias could be reduced by using the FAMz model, instead. 
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Figure 16. Same as fig. ^ for tfie case of FL-rCDM catalogs. 



A quantitative assessment of this bias and its reduction, however, needs to account for the differences between the various 
models and analyses. First of all the linear model used by Nusser et al. (2000) and Branchini et al. (2001) was corrected for 
redshift s-space distortions while our model is not. Second, in this work the v — v comparisons are performed through simple 
linear regressions while the VELMOD analysis of Branchini et al. (2001) and the mode-by-mode comparison of Nusser et al. 
(2000) are less trivial and require additional manipulation of the data that may affect the estimate of /3. While we plan to 
quantify these effects in a future paper we also point out that the tests performed by Branchini et al. (2001) and Nusser et al. 
(2000) have indicated that possible systematic errors on /3 are below the 20 % level. The two previous examples illustrate that 
the results presented in this work are only valid for volume limited and flux limited samples that mimic the PSCz catalog. As 
for any reconstruction method, random and systematic errors need to be estimated again when a new dataset is considered 
or a new model vs. data analysis is implemented. Such error analysis is best done by applying FAM to a set of realistic mock 
catalogs, as we did in this work 

The dependence of the clustering properties on the type of real galaxies implies that most galaxies are biased tracers of 
the underlying mass fluctuations (Loveday et al. 1995, Baker et al. 1998, Hawkins et al. 2001, Norberg et al. 2001). NB have 
shown that a local biasing scheme can be easily incorporated in FAM by assigning to each particle a mass proportional to 
W — {1 + 5)/{l + 5^) where 5 and 5^ are, respectively, the mass and galaxy density contrast in the vicinity of the particle. 
Alternatively, given W one can obtain the mass density field from the galaxy distribution. The mass density field can then 
be sampled by a discrete particle distribution to be used as input in FAM, instead of the original biased particle distribution. 
Any biasing relation is naturally defined in real space; hence an application of FAM on a distribution of galaxies in s-space 
must be done iteratively, if biasing is to be incorporated properly. Nevertheless, several authors have shown that forcing a 
biasing relation in s-space does not introduce large errors (e.g. Verde et al. 1998, Scoccimarro et al. 2001, Szapudi 1998, 
Sigad, Branchini & Dekel 2000, Narayanan et al. 2001). By making reasonable assumptions on the statistical properties of 
the mass density, it is even possible to derive the biasing relation in redshift space directly from the redshift catalog. For 
example, a biasing scheme can be obtained by mapping the probability distribution function (PDF) obtained from the galaxy 
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distribution into a mass PDF which is assumed to be log-normal (Sigad, Branchini & Dekel 2000) or by extensively applying 
hybrid reconstruction techniques (e.g. Narayanan et al. 2001). In future work, we will test FAM using different biasing schemes 
and also using mock galaxy catalogs extracted from "galaxies" identified in A'^-body simulations using semi-analytic models 
(e.g., Kaufi'mann, Nussor & Stoinmctz 1997, Kauffmann et al. 1999, Benson et al. 2000). 

Methods based on the least action principle, like FAM, can serve as a time machine for recovering the particle positions 
and velocities at any early epoch. The particle distribution can then be filtered to obtain the density field as a function of 
time. A backward reconstruction of the large scale structure can be very rewarding. It can provide an estimate for the initial 
density field and this could discriminate among different cosmological scenarios, e.g. by detecting non-Gaussian features in 
the primordial density field or by measuring its power spectrum (Nusser, Dekel, & Yahil 1995, Kolatt et al. 1996, Monaco et 
al. 2000, Narayanan et al. 2001). The performance of FAM as a time machine has only been tested on the scale of the Local 
Group in real space with no biasing (Peebles 1990, 1994), or with a simple halo biasing prescription (Dunn & Lafiamme 1993, 
Branchini & Carlberg 1994). Therefore, there is a need for thorough testing of FAM to establish its reliability at recovering the 
past orbits in redshift space, over scales of cosmological interests and using realistic, time-dependent biasing schemes (Nusser 
& Davis 1994, Fry 1996, Matarrese et al. 1997). While we plan to address these issues in a future paper it is important to note 
that our results, along with those obtained by NB, show that FAMz performs better than linear theory and PIZA (Croft & 
Gaztafiaga 1998) in reconstructing peculiar velocities at the present epoch. It is therefore reasonable to expect that the FAMz 
will also improve over PIZA backward reconstruction and thus will be able to recover the correct initial conditions down to 
scales smaller than ~ 3 h~^Mpc, without enforcing any prior power spectrum like in the Perturbative Least Action procedure 
(Goldberg & Spergel 2000). 

The results of our tests indicate that, once applied to the PSCz catalog, FAMz will able to model the cosmic velocity 
field down to scales of a few Megaparsecs with 1-D random uncertainty of 150 — 200 Kms~^, hence improving over presently 
available models also obtained from the PSCz catalog using either linear theory (Branchini et al. 1999) or the PIZA method 
(Valentine, Saunders & Taylor 2000) . We plan to use this new velocity model to perform v — v comparisons that should return 
an estimate of f2m,o, with an expected la uncertainty of 10 — 15% and, because of the wide range of scales probed by the 
recent data, to provide valuable information on the scale dependence of the biasing relation. 

In this paper, FAM has been applied to systems containing up to iV = 2 • 10* particles. Since the computational cost 

of FAM reconstructions scales as NlogN, the method can also be applied to new generation redshift surveys containing 10^ 
to 10^ objects, like the 2dF and SDSS catalogs. In these new catalogs, however, the surveyed regions span only a fraction of 
the celestial sphere. This constitutes a potential problem for reconstruction methods like FAM which require to be applied to 
spherical catalogs to minimize the dynamical influence of the neglected mass distributed outside the volume of the sample. 
A possible way to proceed is to reconstruct the velocity field over a series of spherical sub-catalogs carved out of the parent 
survey. The sampling within these sub-catalogs will be denser than in the PSCz catalog and the velocity model less affected 
by shot noise. As a result, applying the FAMz procedure to a sot of spherical sub-samples filling the volume of the parent 
catalog would return an unbiased model velocity fields on unprecedented large scales and with 1-D random uncertainties of 
about 100 Kms~^. This would lead to a precise determination of the distribution of mass which could then be cross-correlated 
with absorption features in AGN spectra either due to the interviewing cold, neutral hydrogen (Lya line) or to the warm-hot 
intergalactic medium (Ovi and Oyii lines) to determine the ill-known baryon fraction at z=0. Another possibility would be to 
compare FAM velocity prediction with precise velocity measurements of galaxies in the very local universe, like those obtained 
from the Surface Brightness Fluctuation method (Tonry et al. 2000, 2001). This should allow to measure the value of fim.o 
with an even higher precision. 
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